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Abstract 
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Mathematical models of angiogenesis, pioneered by P. Hahnfeldt, 
are under study. To enrich the dynamics of three models, we intro- 
duced biologically motivated time- varying delays. All models under 
study belong to a special class of nonlinear nonautonomous systems 
with delays. Explicit conditions for the existence of positive global so- 
lutions and the equilibria solutions were obtained. Based on a notion 
of an M-matrix, new results are presented for the global stability of 
the system and were used to prove local stability of one model. For a 
local stability of a second model, the recent result for a Lienard-type 
second-order differential equation with delays was used. It was shown 
that models with delays produce a complex and nontrivial dynamics. 
Some open problems are presented for further studies. 

Keywords Angiogenesis, Nonlinear nonautonomous delay differential 
equations. Global and local stability. Equilibria, M-matrix, Lienard equa- 
tions. MSG 2000: 34K06, 34K20, 34K60. 



1 Introduction 

Angiogenesis, the generation of new blood vessels, is thought to be necessary 
for tumor growth and metastasis [l8]-[20]. This process belongs to a general 
family of tumor-immune interactions. Some recent studies, including mathe- 
matical models of angiogenesis dynamics, are presented in [2], [H], [13], [15]- 
[20], [23] [21], [29], [31] and [32]. 

To incorporate the spatial effects of the diffusion factors that stimulate and 
inhibit angiogenesis, the following two-compartmental model for cancer cells 
and vascular endothelial cells was developed by P. Hahnfeldt in [23] (see, also 

m)- 

dx , , , K(t) 

— = ax{t) m — — 

at x{t) /-|\ 

— = -f^Kit) + Six{t), K{t)) - K{t)) - c{t)Kit), 

where x{t) is the tumor mass and K{t) is a variable carrying capacity, that 
is defined as the effective vascular support provided to the tumor as re- 
flected by the size of the tumor potentially sustainable by it. According to 
P. Hahnfeldt, a stimulator/inhibitor tumor growth dynamics, described by 
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system ([T]), should provide a time dependent carrying capacity under angio- 
genic control and include the distinct mechanisms for angiogenic stimulation 
and inhibition. The dynamics of the second equation is a balance between 
stimulatory and inhibitory effects: first term is the loss of functional vascu- 
lature; the second term corresponds the stimulatory capacity of the tumor. 
The third term reflects endogenous inhibition by either neutralizing endothe- 
lial cell growth factors or inhibition endothelial cell proliferation. These in- 
hibitors are released through the tumor surface (scaling the tumor volume 
to its surface area); thus the major assumption for the Hahnfeldt models is 

I{x,K)/S{x,K) ~ KV, 

where 9 + ^ = 2/3. 

The last term in model ([T]) represents inhibition of tumor vasculature due to 
administered therapy. The following pharmokinetic differential equation 

might be used to model chemo- or radio-therapy mechanisms. Here v{t) is a 
dose given and g is a per capita decay rate of the drug once it is injected. If 
we assume that the drug kills all types of cells |2l], then 



^ = ax{t)\n^-p{t)x{t) 

at x[t) (ij) 

dK ^ ' 

— = ~^iim + S{x{t), Kit)) - Kit)) - cit)Kit). 

at 

To model processes in nature it is frequently required to know system states 
from the past, i.e., models incorporating memory. Depending on the phenom- 
ena under study the after-effects represent duration of some hidden processes, 
for example, time lags of transit through one state to another; transit time 
through compartments, or time lags associated with the growth rates (cell 
division/differentiation time). In any cell growth some cells are inactive and, 
once activated, the cell division is not instantaneous. There are delays in 
cell division. For example, cells can be engaged in an active cell cycle and 
might divide with a fairly regular periodicity. On the other hand, cells can 
be in a resting state (i.e., having dropped out of the cell cycle because of lack 
of positive growth signals) and will experience a delay in resuming the cell 
cycle at the same place that they left. In human cells this delay is usually 
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about 8 hours. One might imagine that in addition to carrying oxygen and 
nutrients to tumor cells in a tumor mass, vascularization may also deliver 
growth hormones (signals) that would stimulate otherwise quiescent cells to 
reenter the cell division cycle and thereby contribute to an increase in tumor 
mass. The Warburg effect on angiogenesis and metastasis is a different mech- 
anism with the after-effects (see, for example, Essentially cancer cells 
continue to use glycolysis much more than regular cells, even when there 
is plenty of oxygen around. This produces an acidic environment (by the 
production of lactic acid from glycolysis) which in turn stimulates angiogen- 
esis. Presumably there is a lag time required to bring the tumor environment 
to the optimal pH for endothelial cell growth. In terms of inhibition, there 
could be time lags between cells receiving an inhibitory signal and actually 
dropping out of the cell cycle. Usually cells will continue through their cycle 
after receiving a signal and all stop at the same spot on the cycle (called a 
restriction point). The inclusion of explicit time lags in the model allows di- 
rect reference to experimentally measurable and/or controllable cell growth 
characteristics (e.g., time required to perform the necessary divisions). In 
general, models with delays produce a complex and nontrivial dynamics: it 
switches stable trajectories into unstable cycles or periodic oscillations. In 
cancer therapy, stability switching is a very important issue in the design of 
a drug protocol (see, for example, [1], [3], [H] and [33]). 

If we assume that the tumor cells enter the mechanisms of angiogenic 
stimulation and inhibition with some delays h{t) < t, then model ([1]) has two 
alternative forms: 
Model 1. 



dx 

Itt 
dK 

~dt 




(3) 



Px{h{t)) - -fK{t) - 6x^^^h{t))K{t) - c{t)K{t). 



Model 2. 



dx 

m 

dK 
It 



ax{t) In— — -p{t)x{t) 
x[t) 

■- l3K{t) - -fK{t) - 6x^/^{h{t))K{t) - c{t)K{t). 




(4) 



Note that a logistic-type model with Richards nonlinearity 
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Model 3. 



dx 

Itt 
dK 

~dt 




(5) 



(3x{h{t)) - -iK{t) - 6x^^^h{t))K{t) - c{t)K{t). 



could also be used for modelling tumor growth dynamics [23]. Here m > 
and m 7^ 1 is a constant that drops an unnatural symmetry of the classical 
logistic curve (m = 1). 

Note that all models without time lags were studied in [13], [IS], [23] and [31]- 
Systems ©-([S]) belong to a wide class of the delay differential equations 



where x G R", A e M"""" and F : R+ x R" x R" ^ R" is a nonlinear and 
continuous vector function. 

Stability analysis of the delay differential equation (E]) is a well-trodden area, 
however, some existing results rely on restrictive conditions, e.g., strict mono- 
tonicity and boundedness of the functions and operators involved, continuity 
of the parameters [12\, [21], [22], [26], [28] and [30]. For example, mod- 
els under study possess non-Lipchitz nonlinearity, thus stabihty analysis for 
models (l3])-(l5l) requires new tools and approaches. Moreover, since the mod- 
els in the paper are based on nonautonomous equations, the application of 
the traditional methods for autonomous equations, such as Laplace trans- 
forms or quasipolynomials, is questionable. 

The purpose of this paper is two-fold: obtain new results in the qualitative 
theory of delay differential equations; apply the results to the analysis of the 
models of angiogenesis dynamics. Firstly, for some general class of nonlinear 
nonautonomous systems with delays we obtained conditions for the existence 
of a global attractor; and based on that result, we proved local stability for 
Model 1. Compared with Model 1, Model 2 has different qualitative features, 
thus for its analysis, we used the theorems recently obtained for a Lienard- 
type second-order differential equation with delays [1]. In the Discussion we 
pose some open problems. 



dx 



A{t)x{t) + F{t,x{t),x{h{t))), t>U 



(6) 
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2 Global Stability Criteria 



We will begin by examining the global stability of a general class of the 
following nonlinear nonautonomous systems with delays. 

n 

Xiit) = -aii{t)xi{t) + '^aij{t)xj{t) + ^fij{t,x{hij{t))), i = l,...,n, (7) 

where \ fij{t,u)\ < bij(t)\u\, aijit) andhijit) are measurable essentially bounded 
functions, fijit,-) is a continuous function, fij{-,u) is a measurable locally 
essentially bounded function and hij{t) < t are measurable functions. Along 
with equations (jTj) we set the initial condition for each to > 

Xi{t) = ipi{t), t<to,i = l,...,n, (8) 

where cpi is a continuous function. 

A solution X(t) = {xi(t), . . . , of problem ©-([H]) is a vector function, 

locally absolutely continuous for t > to, that satisfies almost everywhere 
equation ([7j) on this interval and initial condition ([8]) for t < t^. A unique 
global solution of problem ©-([HD exists. 

Definition. Matrix A is called an M-matrix if aij < 0,i j and one of 
the following equivalent conditions holds: 

-there exists a nonnegative inverse matrix A^^ > 0. 

-the main minors of matrix A are positive numbers. 

Let Aij = supi>ijj \aij{t)\ and Bij = sup^^^^ Ikjit)]. 

Theorem 2.1 Suppose there exist > 0, r > and to > such that 
inf t>tg aii{t) > tti and t — hij{t) < t. The matrix B = {bij} with entries 
bii = Qi — Bii and bij = —Aij — Bij for i ^ j is an M-matrix. Then for any 
solution X(t) = {xi(t), . . . ,Xn(t)}'^ of equation ^ lim(^oo-^(^) = 0. 

Proof. After substitution Xi{t) = e"^^*"*")?/i(t), t > to, where < A < 
miuj ai , system ([7j) has a form 

n 

y,{t) = -{au{t)-X)m{t)+Y,a,,{t)y,{t)+Y,e''~'-'^ 
Hence 



to 
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3=1 



ds. 



Then 



J to 



i(C)-A]dC 



ds 

\a,j{s)\ , , 



\yihijis))\ 



ds 



' aii{s) - A 



^ aii{s) - A 

Set yf = maxto_^<t<6 and = {y\, y^y. Hence 



for to <t <b. We define the matrix C(A) = {cij(A)} with entries 

gAr R . /I . . _^ .Ar i 



Qi(A) 



1 - 



Oi - A 



and Cij{X) 



ai - A 



Hence C{\)Y^ < \X{to)\ for to<t<b.We have hmA^oC(A) = C(0). Since 
matrix B is an M- matrix, then C(0) is also an M- matrix; thus the main 
minors of C(0) are positive. If for some values of parameters tlic elements 
of a matrix are continuous functions then the determinant of this matrix is 
a continuous function. For some A > all main minors of C (A) are positive, 
the latter implies that this matrix is an M-matrix. Assume that parameter is 
fixed A = Aq. Thus for there is an a priory estimate Y'^ < C^^(Ao)|X(to)| 
where the right-hand side does not depend on 6; therefore Y{t) is a bounded 
function. Finally, X{t) = e''^^*~*°^Y{t), then limt-^oo^(^) = 0. The theorem 
is proven. 
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3 Existence of Global Solutions 



In what follows, we assume < t — h{t) < r, x{t) = (pit) for t < to x(to) = xq 
and -ft'(to) = Kq, where p{t) and c{t) are measurable essentially bounded 
functions, h{t) is a measurable function and is a continuous function. 

Definition 3.1 Any solution of problems ^ is a locally absolutely 
continuous function if it satisfies the equation almost everywhere for t > to 
and the initial conditions for t < to ■ 

Theorem 3.1 Suppose that t - h{t) > tq > 0, > 0, x{to) > 

and K{to) > 0. Then systems ^ have unique global solutions 

{x{t), K{t)) positive for t > to- 

Proof. Without loss of generality, we prove the theorem for system ([3]). 
The proof for models 2 and 3 is similar. Suppose for simplicity that to = 0. 
Consider the second equation in system ([3]) for t G [0, To]. This linear equation 
has a form 



Hence K(t) > for t G [0,ro]. Consider the first equation in system (|2]) for 
t G [0, To]. This is a nonlinear ordinary differential equation for the function 
x(t) with a known positive function K{t). Since x(0) > 0, then there exists 
a local solution of this equation. On the interval of the existence of this 
solution we have 



hence on this interval the solution is positive. Suppose that the maximum 
interval of the existence of the solution of this equation is [0,to) for to < tq. 
Since x(t) > then limt^to x(t) = +oo. Therefore there exists < ti < to 
such that 



K{t) = a(t)K{t) + hit), K{0) > 0, 6(t) > 0. 




x(t) > max K(t) for t G [ti,to). 



We have In 



< 0,t G [ti,to). Then 



x(s) 



x{t) < -p{t)x{t) forte [ti,to). 



Finally, 



< x(t) < x(ti)e^*i" W^)!'^' 



for t G [ti,to). 
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It contradicts the assumption that hm^^jp x{t) = +00. Therefore, there ex- 
ists a positive solution of this equation for t G [0, tq]. Similarly, we apply the 
same procedure on the intervals [tq, 2ro], [2ro, Stq], . . . and obtain the global 
positive solution for system (|3]). The theorem is proven. Following the steps 
in this Theorem, it is straightforward to check the existence of the global 
solutions of models (jl]) and (E]). 

In view of Theorem 13.11 we assume that there exists < ti < tq such that 
Ti <t — h(t) < To, e.g., h(t) = t — T. 



4 Models with Constant Rate of Infusion 

Standard chemo- and radio- therapies are typically administered in a con- 
stant dose scheduling, i.e., c(t) = cq and p(t) = po. 



4.1 Stability Analysis for Model 1 



dx 

'dt 
dK 

"dt 



ax{t) In-^^ -Pox{t) 
- /3x{h{t)) - -fK{t) - 5x^'^{h{t))K{t) - coK{t). 




(9) 



Then system 

axln^ — PqX = 

/3x--fK- 6x^/^K -coK = 0, 
has a unique positive equilibrium point 



X* = 

K* = x*e^ 



(10) 



where t] = I3e 



provided that 



/3 > (7 + Co)e 



(11) 
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Let introduce new variables u{t) = In ^ and v{t) = In^^, then system 
has the following exponential form with the trivial equilibrium 



du 
dv 



—au{t) + av{t) 

u{h{t))-v{t) 



rje 



7 - Co + (?7 - 7 - co)e 



2u(h(t)) 

3 



(12) 



At (0, 0) a linearization of system ( 1T2|) has the following form 



du 

~dt 
dv 

'dt 



—au{t) + av{t) 
r/ + 27 + 2co 



u{h{t)) - r]v{t). 



Theorem 4.1 Let t — h{t) < r and condition f|TT]) holds. Then positive 
equilibrium {x*,K*) of system iQ is locally asymptotically stable. 

Proof. To apply Theorem 12.11 we set a matrix 



B 



a 



??+27+2eo 



—a 
V 



Clearly, 



B- 



rj a 

??+27+2co 
3 



a 



2a{ri - 7 - Co) 

Hence B~^ is an M-matrix, provided that inequality f lTT]) holds. By Theorem 
12.11 system f|T2|) is exponentially stable; thus system (|9]) is locally asymptoti- 
cally stable. 



4.2 Stability Analysis for Model 2 

'K{ty 



dx 

H " 
dK 

'dt 



ax{t) In 



x{t) 



Pox{t) 



The system 



{(3 - -i)K(t) - 5x^'^{h{t))K{t) - coKit). 

ax In (-j) — pqx = 

{(3 - -f)K - 6x^/^K - cqK = 0, 



(13) 
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has a unique positive equilibrium point 

x- = ( ^~]~"° )i, ir*=x*e^, (14) 
provided that /3 > 7 + Cq. 

Set M = In ^ and t> = In then system (fT3|) has the following form with 
the trivial equihbrium 

— ^ = —au{t) + av{t) 

at f-N 

Nonhnear system f|T5|) is equivalent to a second order delay differential equa- 
tion 

'^ + a'^ + A(e2"(M*))/3 - 1) = 0; (16) 
dt-^ at 

or, if /i(t) = t , a second order nonlinear ordinary differential equation 

^^„| + ^(,^.../3_i).o. (17) 

where A = a{(3 — 7 — cq). Equations (fT6l)-f|T71) belong to a well-known class 
of a Lienard-type differential equations [6j-[10]. 

^ + /M^ + 9M = 0, (18) 

or 

5 + /H^ + ^7KMt))) = o. (19) 

The following result was obtained by T. Burton ^10] . 

Theorem 4.2 Suppose functions f {u) andg{u) are continuous withf{u) > 
and ug{u) > if u ^ 0. Then the zero solution of equation (1181) is globally 
asymptotically stable if and only if 

r"±00 



^•±00 

/ [f{x) + \g{x)\]dx = ±oo. 
Jo 
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For equation ( 1T7|) we have f{u) = a > 0, g{u) = A(e^"^*''/^ — 1), provided 
A> and u(t) > for all t > 0; thus based on Theorem 14.21 system (fT7|l is 
globally asymptotically stable. 

To prove local stability for Model 2 with delays, we will use the following 
result recently obtained by L. Berezansky et al. ([4j Corollary 5.2). 
To this end, consider the linear equation 

-^ + a{t)- + b{t)u{h{t))=0. (20) 



Lemma 4.1 Suppose a{t) > > 0, t — h(t) < t, and 

m\ 

iim sup — -— < r. 

Then equation [2D^) is exponentially stable. 



Theorem 4.3 Suppose 

t - h{t) <r andO < (3-'j-Co< 1.5r. (21) 
Then the positive equilibrium of system ( 173]) is locally exponentially stable. 

Proof. Stability of the positive equihbrium of system ( IT3ll is equivalent to 
stability of the trivial solution of equation (fT6l) . Linearized equation for (fT6i) 
has the form fl20|) . where a{t) = a and b{t) = |q;(/3 — 7 — cq). Finally, 
application of Lemma 14.11 proves the theorem. 

Remark 4.1 The conditions obtained in Theorem \4.1\ are independent 
of delays, whereas local stability of Model 2 (Theorem \4.3\ ) depends on the 
combination of the magnitude of the delays and some parameters of the model. 

5 Models with Time- Varying Treatment 

Consider general models ([3]) and (jlj), and assume the existence of the limits 

lim p{t) = po and lim c{t) = cq. (22) 

We quote the following lemma [23] . 
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Lemma 5.1 Consider the vector equation 

m 

X{t) = Au{t)X{hu{t)) + (23) 



k=l 



where Ak{t) and F{t) are locally essentially bounded matrix and vector func- 
tions, t — hk{t) < T. If homogeneous equation 

m 

x{t) = Y,Mt)x{hk{t)) 

k=l 

is exponentially stable, and lim^^oo 11-^(^)11 = 0, then for any solution X of 
equation (|23|) 

lim X{t) = 0. 

t—^oo 

Lemma ET] is a simple corollary of the variation of constant formula for solu- 
tions of linear delay differential equations and an exponential estimation for 
the fundamental matrix of this equation. 

Firstly, consider system ([3]) in the form 

^ = ax{t) In (^-^^ - Po^it) + (Po - P{t))x{t) 



dK 

^ = Mhit)) - 7K(t) - 5x^l\h{t))Kit) - coKit) + (co - c(t))K(t). 
dt 



(24) 

After substitution m = In ^ and f = In , system (12^ has a form 

^ = -au{t) + av{t) + - p{t) 

rii, 2u(h(t)) (25) 

^ = ^e«('^W)-''W - 7 - Co - (^ - 7 - Co)e^ + co - c{t), 

where x* and K* are defined by ffTOj) . Linearization for system f l25|) yields 

^ = -au{t) + av{t) + po - p{t) 

^ (26) 

dv r/ + 27 + 2co , , , . , . 

— = ^ u[h[t)) - r/f (t) + Co - c{t). 

Lemma [5.11 implies the following result. 
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Theorem 5.1 Suppose conditions of Theorem 4-i O'^d condition [2^] 
hold. Then the pair of numbers {x*,K*), defined by l[TO\) . is a local attractor 
for the solutions of system 



lim x{t) = X*, lim K{t) = K* . 

t—^oo t—^oo 

Similar calculations yield 

Theorem 5.2 Suppose conditions 021 p and / f^) hold. Then the pair of 
numbers {x*,K*), defined by is a local attractor for the solutions of 

system Q). 

Remark 5.1 Note that the methods applicable for Model 1 yield similar 
results for the logistic-type model ([5]) . 

6 Discussion and Open Problems 

Mathematical modeling and simulation can potentially provide insight into 
the underlying causes of tumor invasion and metastasis, help understand 
clinical observations, and be of use in designing targeted experiments and 
assessing treatment strategies [T], [2] and [R]. Compare to some models 
in literature, where delays were introduced without biological motivation, we 
give a solid justification for the introduction of the delays into the models. All 
models are perturbed either by a constant therapy or by time- varying treat- 
ments. The existence of unique global solutions for models ([3])-([5]) follows 
from Theorem 13.11 For some general class of nonlinear nonautonomous sys- 
tems with delays we obtained conditions for the existence of a global attractor 
(Theorem I2.ip . Based on Theorem 12. H we proved local stability for Model 1. 
For local stability analysis of Model 2 we used recent result for a Lienard-type 
second-order differential equation with delays. Criteria obtained for local at- 
tractivity, are explicit and hence are convenient to apply/verify in practice. 
Note that in Theorems 15. If - 15.21 the study of nonautonomous equations is 
reduced to the study of "limit equations" or asymptotically autonomous sys- 
tems. Note that the analysis of Model 3 could be achieved by using similar 
techniques. 

Finally, we formulate some open problems. 

1. Conjecture Model 1 has a global attractor, provided that a positive 
equilibrium exists. 
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2. Conjecture Model 2 has a global attractor, provided that conditions 
of the Theorem 14.21 are satisfied. 

3. Some additional problems which were not included in the paper: lower 
and upper estimations of the solutions, extinction, existence and asymp- 
totic stability of periodic solutions, oscillation and nonoscillation of the 
solutions about its positive equilibria. 
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